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^ '■ Abstract 

fi \ We illustrate how to apply modern effective field theory techniques and dimensional regularize 
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tion to factorize the various scales that appear in nonrelativistic bound states at finite temperature. 
We focus here on the simplest case: the hydrogen atom. We discuss in detail the interplay of the 
hard, soft and ultrasoft scales of the nonrelativistic system at zero temperature with the additional 



Q-i| scales induced at finite temperature. We also comment on the implications of our results for heavy 



quarkonium bound states in the quark gluon plasma. 
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I. INTRODUCTION 



Finite-temperature effects in atoms were an issue in the early 1980s [l|, 12, y, |j] . The basic 
physics at low temperatures was already understood iu those days Q aud some experiments 

displaying finite-temperature effects were successfully carried out [5|]. The motivation for 
reconsidering this topic is that QED bound states are a good testing ground for heavy 
quarkonium physics |i6j . Indeed quite some number of effective theory techniques, including 
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,112,113|, and 



the use of dimensional regularization, were first tested in QED 7|, la, [ol, 
have now become standard tools in heavy quarkonium physics (see jl^ for a review). The 
behavior of heavy quarkonia states at finite temperature has been believed for a long time to 
be a good probe of the so-called quark gluon plasma |15] (see [l6| for a recent a overview). 



With the advent of current experiments at the Relativistic Heavy Ion Collider (RHIC) and 
the Large Hadron Collider (LHC), precision in the quantification of this phenomenon will 
be necessary, and hence computational tools must be dev elop ed. A number of works in this 
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direction have recently appeared in the literature p/7, 

We present in this paper an efficient way to include finite-temperature effects in nonrel- 
ativistic bound states. We focus here on the simplest of them, namely, the hydrogen atom, 
and make extensive use of nonrelativistic QED (NRQED) 711 and Potential NRQED (pN- 

n n n 

RQED) [8|, |9|, Il0| . Since these effective theories are based on momentum expansions about 
the on-shell condition, which do not exist in Euclidean space, it is compulsory to use the 



real time formalism (see, for instance, 



2^). 



In the hydrogen atom, complications due to hard thermal loops (HTLs) [24, |25|, |26|, |27| 
can be ignored at low temperatures (T -C m, m being the electron mass). This allows one 
to carry out precision calculations in two relevant regimes, namely, when T < E, E being 
the binding energy, and T r^ p » E, where p is the typical momentum of the electron (~ 
inverse Born radius). We critically compare with previous results in the literature. Then we 
move to the high-temperature case T r^ m, which, to our knowledge, has not been studied 
before. We carry out the matching from QED to NRQED at finite temperature and discuss 
the effects of the HTL in the bound state dynamics. 

We distribute the paper as follows. In the next section we review the two effective theories 
mentioned above, which are extremely useful for the description of QED bound states at 
zero temperature, and discuss how they are affected by a finite temperature. In Sections 



Ill, IV and V we address the cases T r^ E, T r^ p and T r^ m respectively. Section VI is 
devoted to a discussion of our results and to some conclusions. Three Appendixes contain 
technical details. 

II. THE HYDROGEN ATOM 

The relevant (energy) scales in the states of principal quantum number n of a hydrogen 
atom at T = are the electron mass m (hard), the inverse Born radius p = ma/n (soft) and 
the binding energy E = — ?7ia^/2n^ (ultrasoft). They satisfy the inequalities m ^ p ^ E, 
which are most conveniently exploited using effective field theories. NRQED is the effective 
theory which exploits the inequality m ^ p,E. It is obtained from QED by integrating out 
momentum scales of order m and is equivalent to it at any desired order in the p/m, E/m 
and a expansions. It reads 



,0 D^ D^ aB |VE| 

2m 8m^ 2m Sm? 

^^^^^ .(D X E - E X D) ^ ^,^^„^ _ I ^ ^^F^^D-F^^ 

D^ = d^ + ieA^ {D^ = d^ — iZeA^, Z is the charge of the nucleus) when acting on ■?/' 
(A^), E (B) is the electric (magnetic) field, and c^, cp, cs and d2 are matchin g co efficients, 



which encode the nonanalytic dependence on the scale m. At 0{a) they read. 
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CZ5 = 1 + - -In- +... 2 

<=. = ! + - + ... (3) 

cs = 1 + - + ... (4) 

TT 
^2 = — + - (5) 

The remaining inequality, p ^ E ys, most conveniently exploited using pNRQED. pN- 
RQED is obtained from NRQED by integrating out energy scales of order p and it is equiv- 
alent to it at any desired order in Ejp and a. It reads 



J 2m Sm-^ 4 

J |Xi — X2I rri'^ o 

Za , Xi — Xo „. X , _i_ , / N 

^ics—^a{-^ \ X V))V^+^(t,xi) 

4m^ |xi — X2|'^ 

The potentials above play the role of matching coefficients, which encode the nonanalytic 

dependence on the scale p. The photon fields in the covariant derivatives contain only energy 

and momentum much smaller than p. This Lagrangian can be written in a manifestly gauge 

invariant form in terms of a wave-function field S'(t,x), which describes an ion of charge 

{Z — l)e and gauge transforms with respect to the center of mass only (it is gauge invariant 



for Z = 1: see 



9| for details). It reads 



^pNRQED — 

+ 



j d^^-^F,,F^^ + j d^^S^t, x) (iDo + ^ + ^ + 

+ f d^xS\t,x)ex-ES{t,x). 

The size of each term above can be obtained using V ~ |x|^^ ~ ma, ido ~ m,a'^ and 
E ~ r/i^a^ (Z ~ 1 will be assumed for the estimates throughout). The leading order terms 
are then in the first line, and produce the well known Coulomb spectrum at 0{mQ^). The 
spectrum at 0{ma^) can easily be calculated from the Lagrangian above, by treating the 
remaining terms as perturbations. The calculation is divided into two parts: (i) a standard 
quantum mechanical calculation of the expectation value of the potentials in the middle line 
between the Coulomb states and (ii) the contribution of the ultrasoft (US) photons, which 
arise from perturbations involving the last term. The former gives. 



(8) 



5'''E,^= ^(f-4d2)K(0)P, 

^"'"^n = c,^ Uj + 1) - Kl + 1) - i) ^nlj\^M3) 



(In) = \nlj)), and the latter 



,-£„. !£f!!fLii + 5_,„2M*'('' 



m b J m? 



(9) 

^|(n|v|r)p(ij;„-ij;.)ln ^ V 

, I -C'n ^r I / 



together with the total width 



rn = $^^«|(Ti|v|r)|2(E„-E,,) . (10) 

where v = —i'V/m and 0n(O) is the wave function at the origin. The correction to the total 
energy is given by 

5En = 5^ En + (5^^E„ , (11) 

in which the /i dependence is canceled between the ultrasoft contribution and the one in cd, 
see (El). 

At finite T, we have to find out how to properly account for the new scale T. The first 
important property, which follows from the Boltzmann distribution, is that fiuctuations of 
energy much larger than T are exponentially suppressed. This implies that for m ^ T the 
same NRQED Lagrangian as for T = can be used: the temperature dependence of the 
hard matching coefficients is exponentially suppressed and hence negligible. It also implies 
that for p ^ T the same pNRQED Lagrangian as for T = can be used: the temperature 
dependence of the potentials is exponentially suppressed and hence negligible. We begin by 
analyzing this case, in which finite-temperature effects are encoded in the ultrasoft photons, 
in the following section. Next we move on to the case m ^ T ^ E. In this case the 
finite-temperature effects must be taken into account in the matching between NRQED and 
pNRQED, and are encoded in temperature-dependent potentials. For T ~ ?7i, the finite- 
temperature effects must already be taken into account in the matching between QED and 
NRQED, and are encoded in the temperature-dependent NRQED matching coefficients and 
in the HTL effective Lagrangian. As in the T = case, we will use the Coulomb gauge for 
calculations in NRQED and pNRQED, and the Feynman gauge for calculations in QED. 



III. THE CASE p > r 

As mentioned before, we can just consider the pNRQED Lagrangian at zero temperature. 
The finite-temperature effects are encoded in the ultrasoft photons, and not in the potentials, 
which remain the same as in the zero-temperature case. Let us count T r^ E and present 
the calculation at order ma^ . If we use the Lagrangian ([6]), there are two contributions to 
the binding energy (and decay width). The first one is given by the photon tadpole arising 
from the kinetic term. It reads {(3 = 1/T), 



(12) 



2' 



3m/3 



The wavy line stands for the tranverse photon propagator (in the Coulomb gauge), and 
the solid line for the nonrelativistic electron propagator. This contribution is bound state 
independent and coincides with the thermal mass shift obtained in direct QED calculations 
[Sj. The second contribution is given by calculating the following ultrasoft loop at finite 
temperature. 



-e^ lim {n\v'Iij{po-Hoy\n) (13) 

pO—fEn 



= — e^ lim y {n\v^\r)Iij{po — Er){r\v^\n). (14) 

r 

The double line indicates that the Coulomb potential is taken into account exactly in the 
propagator, and 

2m |x| J (ztt)'* e'^l'^ol — 1 k^ q — Kq + irj 

We have displayed only the temperature-dependent piece [the temperature-independent 
one has already been given in (^ and (TTUI) ]. If the gauge invariant formulation ([7]) is used 
instead, the whole contribution comes from the last ultrasoft loop. Separating (TT5|) into real 
and imaginary parts, we obtain. 



The intermediate calculations for the imaginary part are presented in Appendix A. We have 
not been able to proceed analytically any further in the general case. We may write down 
our final results for the thermal energy shift and decay width in terms of flT6|) and flT7|) as 



5En = ^^ + e2 5^|(n|v|r)|2$5/,,(E„-E,), (18) 

5r„ = 2e2 5^|(r2|v|r)|2gfJJ,^(E„-E,). (19) 

r 

These final expressions are suitable for numerical treatment. Further analytical results can 
be obtained in the limiting cases E ^ T and E <^T , which we present below. 

A. E^T 

In this case, the real part flT6l) is exponentially suppressed, and hence no temperature- 
dependent contribution to the decay width (iT9l) arises. The imaginary part can be obtained 
by expanding ^ip{iy) for large y in ( TT7I) . 

^^{ty) = ^^{1 + ty) ~ Hy) + j^ + j^ + - (20) 

or alternatively k over q = E^ — Hq in the integrand of (TT5|) . The leading contribution reads 

'' ~ 18/52g' ^^'^ 

so 

Pn = I — Pfi, Pn is the projector onto the subspace of energy En (note that lij{0) = 0). 
This contribution cancels exactly that of the photon tadpole (IT2l) . namely, the first term in 
the right-hand side of ( TT8l) . This cancellation appears to be automatic if one uses the gauge 



invariant Lagrangian (jTj). Either way, the leading nonvanishing contribution conies from the 
third term in fl20l) . 



The matrix element above can be evaluated analytically using the techniques of 28|, |29|]. We 
obtain, 



/n\v- ^^^— -vln) = (nix- — -xln) = —^A(n, I) + l^^B(n, I), (24) 

where 



A{n,l) = - — {F(n,-/-l)-F(-n,-/-l) + 2(n2-/2)[_24(5n2-/2 + l) + 3n(20/ + 15)]}, 

(25) 
and 

B{n,l) = —^—-{F{n,l)-F{-n,l) + 2{n^-{l + lf)[~2A{5n^-{l + lf + l) + 3n{28l + 67)]}, 
bAnmE^ 

(26) 
with 

F{n,l) = {n + l + 2){n + l + l)[^^^±^±^^^^^—^ + A{2n-in (27) 

The details of this computation are given in Appendix (A. 2). Note the strong dependence of 
the expression above on the principal quantum number ~ n^ /ifjv'a^. Let us then summarize 
our final results for the thermal energy shift and decay width in this case as 



47r^a , , P„ , , A. .^ ( , r? 



5r„ = 0. (29) 



B. E <SCT 

In this case, the real part can be easily evaluated by expanding the exponential. At leading 
order in this expansion, it leads to an additional temperature-dependent decay width for all 
the states. 
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''" - w- (»' 

The total width is obtained by summing the T = contribution ( fTOl) to the expression 
above. The imaginary part is obtained by doing the y — > expansion in (ITTj) . 

^ij{zy) = -j + 0{y^). (31) 

Ahernatively, one may expand q = En ~ Hq over k in flTSl) . Then the Bethe logarithms from 
(fT7|l cancel out against those of the zero temperature contribution Qj, and we get for the 
whole ultrasoft contribution, 

3 ivr 6 m"' 

The total binding energy is obtained from (fTTj) using the expression above for 6^^ En and 
(IH]) for 6^ En- Alternatively, we may summarize our final results for the thermal energy shift 
and decay width in this case as 

c-^ ctTT 2« v-^ , , , , ,,r,,^ ^ , ,, , 27r , ,, ,^ f ,l3'ma.o\. , , 

3/3^2 \ \ n- 



*r„ = ;;^(i+o(^)). (34) 



IV. m:^T:^ E 

In this case finite-temperature effects are expected to modify the potential, which might 
in principle give rise to qualitatively different effects. However, for QED at energies below 
the electron mass the vacuum polarization effects are suppressed by even powers of m, and 
hence the full Aq propagator in the Coulomb gauge is not sensitive to the temperature 
(up to high orders in T/m (~ T^/m^)). Finite-temperature effects enter only through 
the tranverse photon propagators. Since the coupling of these photons to nonrelativistic 
electrons is suppressed by powers of l/m, the finite-temperature effects modify only the 
1/m corrections and, hence, the Coulomb potential remains as the leading order term. This 
implies that the gross features of the hydrogen atom spectrum will be kept the same for 
temperatures smaller than the electron mass. We proceed then to the matching between 



NRQED and pNRQED at finite T. At T = the matching is trivial in the electron sector, 
since this sector is insensitive to the momentum transfer (to transfer momentum one needs 
the nucleus), the soft scale to be integrated out. At T 7^ the temperature is a scale to be 
integrated out and the matching becomes nontrivial in this sector. If we count T ~ p, for a 
calculation at order ma^ we need the contributions of the following diagrams: 



.^— ,[ln- + 7-ln2 + -](po--)-^^^, (35) 



macp 



p^ 



3^2 



{Po-TT- 



tan c 



2m 30m^/9 



4' 



(36) 






(37) 



ma 
3m/32' 



(38) 



relativistic 






(39) 



The last diagram comes from the ip'^Ty^ip/Sm^ term in the Lagrangian, which contains a 
piece with two derivatives and two A fields. Other possible diagrams either are of higher 
order or give zero. 

The first diagram has an infrared (IR) divergence. We have followed the same prescrip- 
tions as in the T = case. We have regulated it in dimensional regularization (DR) and 



10 



used the modified minimal subtraction scheme [MS). When one will eventually make cal- 
culations in pNRQED one must regulate the ultraviolet (UV) divergences which will appear 
there in DR and use the same subtraction scheme. The subtraction point dependence will 
then cancel out in all observables and the finite pieces will be consistently calculated (see 
|9|] for detailed discussions in the T = case). 

The matching in the electron-nucleus sector (i.e., the calculation of the potentials) reduces 
to the calculation of the following vertex diagrams: 



where 



^— ^[ln- + 7-ln2 + -], (40) 



^ = ^ ^, (41) 





Cp ^^^-^ cp 



6m2/52- ^ ^ 



The dashed line stands for the Aq photon propagator (in the Coulomb gauge). As before, 
other possible diagrams are zero or of higher order. Putting all these together we obtain the 
following temperature dependent corrections to the pNRQED Lagrangian: 

SL{T) = I rf3x[^[ln(/?/i) - ln2 + ^ - ln(27r) - 7][^|^ + A^+5oV^] - 
TTO; _ , , _ , 2TTac% ^'V'^p'^'V^p . ,, „ ,, , an , , , ,, 



+ f ^^xirf^X2iV+iV(t,X2)( ^ fr^?^""^ | ^^^(^,xO + (43) 
J bm'^p'^\x.i — X2I 

+ |A | ^" , VV^+VV^(t,Xi)[ln(/?/i) - ln2 + ^ - ln(27r) - 7]), 

O TlTfl X]^ — ^2 ^ 

which can be cast into a much simpler form by using the following field redefinition: 
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2a r, /3/i T O ^1 A '^^^ 



*-(' + w''"i?+^-'"^-^6)^-(w^>^'- (^^) 



^^(^) = /''''=|-6£^^*'^*'+(35^W'"*l + 



7- (ml — , , , , 



+ / ci3x,rf3x2Ar+Ar(t,X2)(-^;3^(ln(^)+7-ln2 + -)53(xi-X2))^+V'(i,Xi). (45) 



In order to calculate the spectrum at the desired order we only have to sandwich the poten- 
tials between the states and calculate the US contribution (and, of course, take into account 
the relevant mass shifts in ( l45l) ). The first contribution gives 

The US contribution corresponds exactly to the diagram (iT3l) , but it has to be calculated 
taking into account that it contains now only energies much smaller than T. In this case the 
Boltzmann factor can be expanded. This may (and will) introduce UV divergences, which 



as mentioned before, must be regulated in DR and MS subtracted in order to be consistent 
with the calculation of the potential. We obtain 



^UM) = (l7r^S,,i\nf- + - - In 2) + 0(g^/3^), (47) 

ovr^ |g| 

SJjj(g) gives a contribution to the binding energy which exactly cancels that of the T = 
piece OH]). Then the total binding energy is obtained by adding to fHB]) the T = soft 
contribution ([H]). ^Iij{q) gives a contribution to the decay width which coincides with fl30p 
at leading order in the ma/T expansion. This contribution is parametrically larger than the 
zero-temperature decay width fITOl) . Notice also that in the limit p ^ T the binding energy 
( 146|) reduces to (l32l) . as it should. We may summarize our final results for the thermal energy 
shift and decay width in this case as 



^^ cc^ 7i"o^ 2q; \— \ , , , , ,,9,^ „ w, / 27r , , ,, , 

3/3n2 ^T 



AZ'^a^ 2a v^ , , , , ^ ,n, ^ ^ , ,, , 



The results above are accurate up to corrections of order ma^ for temperatures T ~ ma/n. 
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V. THE CASE mr^T 

For temperatures of the order of the electron mass, electron-positron pairs are created in 
the thermal bath, which are expected to destabilize the hydrogen atom. In order to make 
this expectation quantitative, we will integrate out the scale m ^ T. In the photon sector, 
this will induce a mass dependent HTL effective Lagrangian. In the electron sector, not 
only will the NRQED matching coefficients now depend on T, but also new nonlocal terms 
appear. Let us analyze these two sectors in the following. 

A. Matching QED to NRQED+HTL 

1. The photon sector 

The HTL effective Lagrangian will be obtained from the vacuum polarization, by standard 
techniques 22|. Rather than depending on the single scale eT, as in the massless case, the 
HTL effective Lagrangian is now expected to have a nontrivial dependence on m/5. In fact, 
this brings in a new qualitative feature: the angular integration appearing in the massless 



case becomes a full three-parameter integration 3l|]. In order to illustrate it, let us focus on 
the longitudinal component of the retarded self-energy, which will be the only one needed 
later on. Using the fact that po, p <^ m, T and expanding them accordingly we arrive at 

H?°(p) = (-z)2e2 / —7^ "^^^^ o- (50) 

Y'^ v^i?T^ ^ "I J 
Note that when m = the integral over fc = |k| decouples from the angular integration 
and can be carried out analytically. For m ^ 0, however, the integral over k remains in the 



effective theory. If we write it in terms of w := k/yk^+m? {w G [0, 1)), it is clear that the 
HTL effective Lagrangian for the photons can be obtained from the one in the massless case 



(see, for instance. 



22| ) by doing the following substitutions: 

k — ^ w (51) 
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(1-^2)2 I eVI^^ + 1 
k = k/A; and dVt is the integration measure of the sohd angle. 

2. The electron sector (NRQED) 

We have just seen that the photon sector at finite temperature is quahtatively different 
from the zero-temperature one. Indeed, in the former case a nonlocal HTL effective La- 
grangian is produced, which is much more important than the l/irP' suppressed terms that 
arise at zero temperature (last term in (1)). The question is then whether in the electron 
sector something similar will also happen. In order to find out, we match QED to NRQED 
in this sector as follows. We calculate the two-point Green function of an electron with 
momentum p^ and sandwich it between P^ = ^-y^. Then we make the change Pq = m + ko, 
p = k and expand for ^o ~ ^ and k small. We will find that, unlike the photon sector, 
the expansion is local. Then it will be possible to identify 6Z^{k), the matching coefficient 



of the nonrelativistic field (P+^ = \/ Z^{Q)ip + 0{l/m), where \E' stands for the relativistic 
Dirac spinor field of the electron), and Q{k), the NRQED self-energy. 



iSZ^k) ^ iQ{k) ^___ 




ko - k'^/2m {ko - P/2m)2 



In the real time formalism, the propagators consist of a sum of the zero temperature part 
and the thermal part, which will be proportional to ns for photons and np for electrons 
{nB(F) are the Bose (Fermi) Boltzmann distributions, ub^f) = l/(e'^''^°' =F 1)). If we consider 
just the contribution of hb (np), we are taking into account the thermal fiuctuations of the 
photons (electrons). It is important to note that in the diagrams we will consider it will never 
appear UbUp terms because of kinematic constraints (we will never have an internal electron 
on shell and an internal photon on shell). Hence we can write 6Z^{k) = 5Z^{k) + 6Z^{k) 
and Q{k) = Q^{k) + Q^{k). 

Let us first consider the contributions from the thermal fiuctuations of the photon to the 
electron self-energy. We obtain 
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e^(A;) = -I^-I^ + of — 

6Z^ik) = —ll^+ \- +0( — ). 53 

Note that Q^{k) corresponds to a thermal mass shift 6^m = vra/Sm/?^ for the electron. 
6Zl^ is IR divergent because of I a'- 

0D_i is the solid angle in D — 1 space dimensions, D = 4 + e, e ^ 0, and we have used the 



MS subtraction scheme. 6Z^ will be relevant for the calculation presented in Appendix C 






For the thermal fluctuations of the electrons we find a similar result. Q (k) gives rise to 
the following thermal mass shift 



S^m = ^himP) - ^^3^^, (55) 



p Aam 2ag{m(3) 

^ ' h{mp) 

TT vrmp 

h{mf3) and g{mf3) are defined in Appendix B. Note that 5^m above goes to zero ex- 
ponentially ii m ^ T. 5ZK{k) is simply related to the thermal mass shift 5ZK{k) = 
-6^m/m + 0{k^/m^). 

In principle we should have taken into account the doubling of degrees of freedoms in this 
calculation, as we did in the photon sector. However, the off-diagonal components of the 
self-energy vanish for the same kinematical reasons that forbid terms proportional to hbUp 
above. Hence, the self-energy is diagonal and we can safely ignore the doubling. 

In view of the above results, we may wonder if any QED Green function involving electrons 
will match to local NRQED operators, as is the case of the two-point function, or new 
nonlocal HTL vertices will arise. Let us then analyze the vertex (three-point function with 
two-electron and one-photon legs) next. The calculations are presented in detail in Appendix 
C, here we summarize only the more important results. If we just consider the thermal 
fluctuations of the photons, the vertex can indeed be matched to local NRQED operators. 
In Appendix C we display the modifications of cd and cs in ([2]) due to temperature in the 
case T ^ m as an example. However, if we take into account the thermal fluctuations 
of the electrons we get a nonlocal vertex (see (]C8I) in Appendix C). This vertex is of the 
same size as the tree level contribution when the momentum transfer q ~ ma (the typical 
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moinentum transfer of the bound state at zero temperature), and it is suppressed only by a 
factor e when q ~ me (the scale of the Debye mass). Hence, it turns out to be much more 
important than the local contributions arising from the thermal fluctuations of the photons 
when T ~ -m. Nevertheless, it goes exponentially to zero when T <^m. 

The effective theory for a nonrelativistic electron in a thermal bath of T ~ m 
(NRQED+HTL) lies then in an intermediate situation between the case T <^ m (NRQED), 
in which all contributions are local, and the massless case, in which all contributions are 
nonlocal (HTL). 

B. Matching NRQED to pNRQED with HTLs 

We shall restrict ourselves to the leading order contributions. The matching is then 
analogous to the T = case, which leads to the Coulomb potential, but using the HTL 
propagator for the (Aq) photons. The latter can be obtained from the retarded self-energy 



501) by a standard procedure (see, for instance, |32|]). It reads 



^ / X _ „■/ 1 iWagjmp) 

where we have used Po *^ P? P = |p|) and r7i|, and g{mj3) are defined in Appendix B. By 
Fourier transforming, we obtain the following real space potential: 

where 

(58) 



:^ dzz 
0(x) = 2 



{z^ + l? 



s\.mzx] 



zx 
Unlike for the T = case, now the Aq photon propagates over arbitrary times, which, 

together with the fact that its propagator contains scales, implies that contributions to the 
self-energies of both the electron and the nucleus arise. These read, 

5m = -^^ - ,!^!^M) (59) 

for the electron, and the same expression multiplied by Z^ for the nucleus. In order to per- 
form this calculation we need, in principle, Aii(p, po) for any kinematical region. However, 



due to the fact that Aii{p,po) = Aii{p, —po) (see, for instance, |22|]) we have 
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/ 



dpo -—Aii{p,po) = 7rAii(p,0) , (60) 

Po + iV 



and hence the expression (156|) is enough to carry out the calculation. Formulas (!58l) and (159|) 
are analogous to the results obtained in [17] for QCD, which we recover in the m ^ limit by 
setting Z = 1, and changing e^ -^ g"^ and the group factors in m|,, namely, 1 —^ Ca + Nf/2. 
Notice, however, that our calculation is much simpler: only one tree-level and one one-loop 



diagram need to be calculated, instead of the five one-loop diagrams needed in ref. [l7|. It is 
important to realize that (158|) has an imaginary part. The HTLs induce the scale mj:, ~ eT, 
which for m ^ T dominates over the typical momentum scale of the bound state at T = 0, 
p ~ ma, and hence dramatic changes in the bound-state dynamics are expected to occur. 
Indeed, if p ~ rriD then the imaginary part of the potential is more important than the 
real one and no bound state is expected to survive. The typical momentum for which the 
imaginary part becomes of the same order as the real one is p ~ (16aY^^ g^^^ {mP)T =: m^. 
rrid may be considered as a new dynamical scale in the system, which is parametrically 
larger than mo- Notice that both m^ and rra^ have a nonanalytic behavior in m/T: when 
T becomes smaller than m they go exponentially to zero. The leading behavior of rriD and 
rrid for T < m reads 



ml ~ SaJ^e-'"^ (61) 



2/3 



/lQam\ ^ 



2m/3 



Note that m^ is exponentially larger than tti^j. This allows us to get more explicit expressions 
for the energy shift and the decay width of the bound state in the case m^ "^ f^d ^ P ~ 1/r, 
upon expanding ( l58ll on mo, and using the asymptotic expressions for tjid and rrid above. 



arriD^Z — 1)^ Zam 



SEr^i = -^^^^^- '- - ^^^^{nl\r\nl), (63) 



6T^i = 2{Z-lfaJ—^ ^(n/|r2(lnm^r + 7-4/3)|nO. (64) 

The expressions above hold up to corrections 0{{'mj:,n'^/may). {nl\r\nl) = [Sn^ — l{l + 
iy\/2Zma and (n/|r^|n/) = [5n^ + 1 — 3/(/ + ly^n^ /2{Zmaf' can be found in standard 
textbooks, and 
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(n/|r^lnr|n/) 



9n 
2{Zmay 

In. " 



(n-l-l)\ 



n-l~l 



r(2/ + 5 + r) 



2Zma 



(n + l)\ ^^ r(r + l)r2(n-/-r)r2(5 + /-n + r) 

^ ^ r=n—l—4: V / \ / \ / 

+ V^(2/ + 5 + r) + 2^(4) - 2^(5 + I - n + r)). 



(65) 



which may be obtained using, for instance, the techniques of ref. 33|. 

For Z 7^ 1 it is interesting to observe that the system develops a decay width that is 
not exponentially suppressed (first term in fl6^ ). This is because a charged ion will tend to 
capture electrons from the thermal bath to decay into a less charged ion and eventually into 
a multielectronic atom. Let us focus in the Z = 1 case. For n large enough, namely, when 
ma/n ~ mD, the approximation that leads to (IMl) above fails. However, much before, when 
ma/n ~ md, the states n will melt, namely, their decay width will become of the same order 
as the binding energy. Therefore the expressions in flMl) are appropriated for T < m as far 
as it makes sense to speak about states with a narrow width. 



n 


Td (keV) 


iriD (keV) 


md (keV) 


1 


60.4 


0.703 


3.73 


2 


50.1 


0.284 


1.86 


3 


45.6 


0.167 


1.24 


4 


42.9 


0.114 


0.932 


5 


40.9 


0.0842 


0.746 



TABLE I: Dissociation temperature, Debye mass m£), and dissociation scale tti^, as a function 
of the principal quantum number n. The dissociation (melting) temperature is defined as the 
temperature for which the dissociation scale m^ equals the soft scale mZa/n. Note that rriD is 
smaller than the soft scale but much bigger than the ultrasoft scale m{Za)'^ /n'^, which is consistent 
with our assumptions. 



VI. DISCUSSION AND CONCLUSIONS 

We have developed a formalism which allows us to efficiently factorize the various scales 
appearing in nonrelativistic bound states at finite temperature. It makes use of dimensional 
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regularization and of the known Effective Field Theories both for nonrelativistic bound states 
(NRQED,pNRQED) and for finite-temperature systems (HTLs). We have focused on the 
hydrogen atom. 

For T ^ m we have calculated the finite-temperature effects to the binding energy and 
;he decay width to a precision equivalent to O(ma^). We agree with the early results of 

1 nnn 

fj, but disagree with others [2|, y, |J]. It is interesting to recall how the finite-temperature 
effects were experimentally observed in atoms in the early 1980's 5|. Since En ~ ma'^/n'^, 
even if for the ground state Ei ^ T, there will always be n's, n ^ 1, for which En ^ T. 
For the ground state, finite-temperature effects may be very small (given by ( l23l) and ( 12^ ) 
but for highly excited states the thermal mass shift flT^ must arise. Then transitions from 
highly excited states to the ground state are sensitive to the thermal mass shift. 

For T ~ m we have restricted ourselves to discussing the dominant effects due to finite- 
temperature. In the photon sector, we have described how to obtain the HTL effective 
Lagrangian for a finite electron mass. It requires the introduction of an extra integral 
in addition to the solid angle one. In the electron sector we have seen that in addition 
to temperature-dependent NRQED matching coefficients, new nonlocal (HTL-like) terms 
arise. We have calculated the potential at leading order, which develops an imaginary part. 
The massless limit of this potential agrees with the Abelian limit of the one obtained in 



171]. The imaginary part dominates over the real one for momentum transfer smaller than 
iTLd ~ e^'^Tg^'^(mP). For T < m, g{mP) increases exponentially from zero when T increases. 
Then, for a given bound state, there will always be a temperature for which the soft scale 
equals rrid, and hence the imaginary part (decay width) equals the real part (energy). We 
call this temperature dissociation temperature and have calculated it in Table 1 for the 
lower-lying states. For temperatures higher than the dissociation temperature, it does not 
make much sense to speak about a bound state any longer. 

We then get the following picture of a hydrogen atom in the ground state when heated 
from T = to T ~ -m. The effects are very small until T ~ ma'^. Then it starts developing a 
width ~ Ta^, which increases with temperature but remains much smaller than the binding 
energy until T r^ m. Then, the width starts increasing exponentially and the hydrogen atom 
ceases to exist. 

From our results we can infer some qualitative features of heavy quarkonium systems 
in the weak coupling regime (i.e., when the binding is due to a Coulomb-type potential) 
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at finite temperature. These states satisfy Aqcd ^ ma'^, Aqcd being a typical hadronic 
scale. Aqcd affects at most the next-to-leading-order corrections, and hence these states 
are expected to be rather insensitive to the QCD deconfinement phase transition. When 
the temperature overcomes the ultrasoft scale (T > ma"^), a decay width proportional to 
the temperature will be developed, analogously to the hydrogen atom. As the temperature 
increases further, gluons and light quarks will induce a HTL imaginary part in the potential 



17l |. which will become comparable to the real part when T ~ mas ■ No bound state is 
expected to survive beyond that temperature. One should keep in mind, however, that only 
the ground states of bottomonium (T(IS') and rib), and to a lesser extent of charmonium 



{J/4' and ?7c), are likely to be in the weak coupling regime [34 1. 
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APPENDIX A: CALCULATIONS IN PNRQED 

1. Self-energy 

We proceed to the detailed calculation of the self-energy in pNRQED. It is convenient to 
separate it into real and imaginary parts. The real part is immediate to obtain and has been 
given in (iT6ll . so we will focus on the imaginary part. We expand the Boltzmann distribution 
function in flT^ as follows 



o(3k _i kA^nd(3 
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lE^j!?^-"^'. (Ai) 



and get 



oo 



S/., = -Ij^ ^(e'^^l'^li^iMlgl) - e--P\'^\E*{nP\q\)), (A2) 

where -E'i(x) = J dte~^/t and E*{x) = —P J ,dte~^/t {P stands for the principal value). 
Now we use the following property of the above functions |30| , 

r ^?T7F ^^ = ^ [e"^i («^) - e~"^* («^)] ' (A3) 

Jq ci + 1 Z 

and get 



V(e"'^l^lEi(n/5|g|)-e-"^l''li?*(n/5|g|))=2 / dttcos(/?|g|t) V ^— ^ . (A4) 

n n 

The sum can be carried out using complex variable techniques. We obtain 

_ 2 ^^ /-"^ cos(/J|g|t)rft 
^^^^ " 3 (27r)2 X ttanh(vrt) ^^^^^^^^^ ^^^ (^5) 



Finally, the integral yields, 



2 5ijq 27r «> //^/^kkx /.^n 



where '?/'(x) = r'(a;)/r(a;) 



2. Computation of (124]) 

We derive here the result displayed in fl24p . A similar computation has been done in the 
past for QCD 28(]. However, (12^ cannot be obtained by just taking the Abelian limit of 
the QCD result. The latter is singular because it does not contain the projector P„. We 
will proceed in a way analogous to [28j] , but keeping E ^ En'vn the terms in which the limit 
E ^ En does not exist. We set E = E^ + \ with A — > 0. If we drop Pn in the numerator of 
([MD we get 

{nyjj^n\n) = ^-^ + g{n, I) + 0(A). (A7) 

If we included Pn in the numerator we get an exact cancellation of /(n, /), so that the limit 
A ^ can be taken safely: 
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1 f 1 

{ny ri\n) = / d^x £y{ny\x){x\ \y){y\ri\n), (A8) 



we use r^\x) = x^\x) and the following formula 281] 



1 °° j 

z=o ^ 



with 



(AlO) 

Note that ii E = En then Zma/k = n, so we will have a pole at s = n — / — 1. For the 
angular integration we take into account that x^yi = xy{^^) and combine it with P; using 



.PA.) = (' + iW..M + 'f.-.M. (Au) 



We get finally 



^' H-E'^' n{2kof{2l + l){n + l)\j^^{s + l-Zma/k){s + 2l-l)\ ^ ' 
2mk{l + l){n - I - l)\ ^ (s!(/(/ + l,s))2 



n(2A;o)^(2/ + l)(n + /)! ^ (s + / + 2 - Zma/A;)(s + 21 + ?,)\' 



with 



J(/i, s) = 2A;o / rfx(2A;ox)'+^(2A;x)^-('=+'=°)"L^'+/_i(2fcoa;)i:f +'(2A;x). (A13) 

We define fco so that, En = —kQ/2m. For terms in the sum that are well defined when k ^ k^ 
we can just put k = k^. For the terms in the sum that are singular we have to expand for 
small A and then subtract the singular part as we mentioned before. This is indeed what 
the introduction of P^ = 1 — Pn does. In order to demonstrate it let us look at the Pn part. 

{n\r' ^ ri\n) = "t / d^xd^y{ny\x){x\Pn\y){y\ri\n). (A14) 

We will proceed in an analogous way as in the calculation with no projector above, so that 
the cancellation will become apparent: 
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ra-1 ,• 



(x\Pn\y) = Y.{2l + l)Gi{x,y)Pi{^, (A15) 

1=0 ^y 



where 



TT mn (n + /)! 

is similar to (^^(x, y; —k"^ /2m) above , but the summation for s is restricted to singular terms. 
By comparing ( JA16I) and (1A10I) we can easily see that the 1/A terms cancel even before doing 
the radial integration. Since the P„ part is proportional to 1/A (no finite pieces), we only 
have to calculate the finite contribution (in an expansion in A) to the part with no projector. 
This can be easily obtained by expanding k about k^ in (JAISP (recall that the derivative of 
a Laguerre polynomial is a Laguerre polynomial). The computation can be terminated in 
an analytic form using, 



Jo 






r=0 



n ~ r / \n — r / \ r 



APPENDIX B: INTEGRALS IN TERMS OF SPECIAL FUNCTIONS 

We give here the definitions of various functions appearing in the paper in terms of one 
parameter integrals and provide expressions in terms of special functions. The Debye mass 
can be expressed as 

ml := ^e\2f{mP) + h{mP)), (Bl) 

where 

/(m/3) := i, / dk ^ ^' , = - f (-l)"^i^, (B2) 

and g{rnP) as 
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g{mP) := ^^ 



dk- 



,l3Vk 



^L^^ = m/51n(l + e-^) + U,{-e-^) + ^ 



2/32 



m^(3 



IQ gAJVt--t-'"- _|_ I o 2 

Kq{x) and -ft'i(a;) are Bessel functions and Li2{x) the dilogarithmic function. 



(B4) 



APPENDIX C: MATCHING THE VERTEX FUNCTION 

1. Matching QED to NRQED+HTL 

In order to carry out the matching for the vertex function (F) we have to deal with 
the doubhng of degrees of freedom. There are three external particles in the vertex, and 
each one can be of type 1 or type 2 (following the notation of [22']), so F is a tensor with 
eight components. But, because of kinematic constraints it cannot happen that there is 
and internal photon on shell and an internal electron on shell at the same time, so the only 
components that are nonvanishing are 111, 121, 212 and 222 (the middle index corresponds 
to the photon). If we take into account that the matrix elements of the propagators in the 
real-time formalism are not independent, we obtain Fm = F222 and F121 = F2i2- Notice also 
that, for the physics of an atom with an infinitely heavy nucleus, the only components that 
have a contribution at first order are 111 and 212. 

As we did with the self-energy, we calculate first the contribution from the thermal 
photons. In this case F212 = 0, so we only have to calculate the 111 component. The 
calculation is done by matching three-point Green functions in QED and NRQED: 




(CI) 




(C2) 



The first row represents QED diagrams (all of them are sandwiched between the projectors 
^y^), and the second one represents NRQED diagrams. We find 
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5e = 0, (C3) 

CD = l + ^(ln^ + 7-ln2 + ^ + 0(T2)), (C4) 

6cs = + O{aT^). (C5) 

The finite-temperature contribution to Cp is irrelevant at first order for a hydrogen atom since 
the corresponding operator contains tranverse photons only. We have restricted ourselves 
to calculate the leading order contribution to cd and cs in the limit of T -C m, which is 
enough for illustration purposes. This is also justified because for T ~ m the neglected 
contributions, as well as the ones taken into account, produce modifications of the spectrum 
of order ma^, whereas we will see in the following that there are contributions from the 
thermal fiuctuations of the electrons at order ma^, which will modify the physics of the 
hydrogen atom in a much more profound way. 

We focus next in the contributions from the thermal electrons. The leading order contri- 



bution in QED comes from the second diagram in fICip 



1 





dw r dQ 1 



^Fiii = iSir'^ame / ^^ / 7^(1 - — ^ )S{Qo - qw), (C6) 



dw f dn 



6r2i2 = i4:7c'^ame I ^„ / 7^^^ ^ )6{qo-qw), (C7) 



(w = k/y/k^ + m^, k being the momentum circulating in the loop). Note that this contribu- 
tion is nonlocal and cannot be matched to any of the NRQED operators. It can be matched 
to the following nonlocal operator 

SL= f d^wf{w)wE^6{ido - zwV)^+V^, (C8) 



where 



/H = ^ — -13;^ 2n iw?U~^ • (^9) 

(E is the electric field). This operator becomes as important as the leading order Lagrangian 
when q ~ ma, and it is suppressed only by e when q ~ me, the scale of the Debye mass. 
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Hence the thermal fluctuations of the electrons have a bigger impact in the NRQED La- 
grangian than any of the relativistic or radiative corrections. 

2. Matching to pNRQED and cancellation of the scale dependence 

The nonlocal vertex above can easily be matched to pNRQED at tree level by expanding 
the energy over the three-momentum. At leading order we have 

SLpNRQED = I rf'w/(u;)wE^5(-zwV)^+^. (CIO) 

J wV 

This vertex is IR divergent, so in the calculations in pNRQED there should appear an 
ultraviolet divergence, in order to get a cancellation of the /i dependency. It indeed appears 
in a diagram of the type. 




(Cll) 



where the internal lines are now nonrelativistic propagators for the electrons and HTL 
propagators for the photons^. 

For simplicity, let us check this cancellation in a specific piece of the tensor vertex (for 
the remaining pieces it will be analogous). We focus on ^Fm, in the case go -^ 0, and take 
into account only the temperature-dependent part in one of the electron propagators and 
the zero-temperature part in the other one, which will be enough for illustration purposes. 
Let us call it F*. 

From the NRQED matching we have (from the first term in (lC6p . by taking go = and 



undoing the change of variable w = k/ yW^+m?) 



(C12) 



Since we are interested in only the IR divergent behavior, we may substitute the integrand 
by the following regulated expression 



^ Note that Coulomb resummations can be ignored at the scale of the Debye mass (me), since they only 
become important for momentum transfer of the order ma or smaller. 
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r 7 — ,,, ,, ^ r/i / ^ = -7 — TT^-^ -+n— . C13 

Any calculation in pNRQED+HTL involving the contribution above will also involve the 
diagram fICllI) with nonrelativistic propagators for the electrons and HTL propagators for 
the (Aq) photons. Let us take into account only the temperature-dependent part in one of 
the two electron propagators and the zero temperature part in the other one, in accordance 
with the evaluation of T* above, and call the corresponding contribution T*. We have. 



d'k , {-27r)6{ko-£)_ i 

(-^g) .^i^^^.i , r (-^e) : un:iA2 -(-'^e)An (|p - k\,po - k^) 



(C14) 

Due to the 5 function and to the fact that po ~ p^/2m ^ p, we can use the expression f l56|) 
for All. We focus on the UV behavior of the expression above, since we are only interested 
in identifying the /i dependence, which should cancel that of fIClSp . We can then neglect 
the imaginary part of An, which leads to finite expressions, and approximate 



T.* 3 ' ^'k 1 1 



(27r)3 e^™ + 1 -(qk) + irj (k)^ + m]^ 

Tie'm 1 r dkk^+' 

(C15) 



"D 



3^™ + l)|q| (27r)2 7o P + m' 
' -- + -ln ' 



(e^™ + l)|q| (27r)2 V e 2 \ml 

In the second relation we have carried out the angular integration and introduced DR (ne- 
glecting e in the finite pieces). If we add T* to T* we see that the fi dependence indeed 
cancels, as it should. 
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